Evidence for Macroscopic Quantum Tunneling of Phase Slips in Long 
One-Dimensional Superconducting Al Wires 
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Quantum phase slips have received much attention due to their relevance to superfluids in reduced 
dimensions and to models of cosmic string production in the Early Universe. Their establishment in 
one-dimensional superconductors has remained controversial. Here we study the nonlinear voltage- 
current characteristics and linear resistance in long superconducting Al wires with lateral dimensions 
~ 5 nm.We find that, in a magnetic field and at temperatures well below the superconducting 
transition, the observed behaviors can be described by the non-classical, macroscopic quantum 
tunneling of phase slips, and are inconsistent with the thermal-activation of phase slips. 

PACS numbers: 74.78.Na, 74.25.Fy, 74.25. Ha, 74.40.+k 



What is the resistance of a superconductor below its 
superconducting transition temperature? For a three- 
dimensional superconductor, the answer is the obvious 
one - zero. But in one-dimension (ID), by which we 
mean very thin very long wires, quantum fluctuations 
destroy the zero resistance state. Phase slips - small su- 
perconducting regions that become normal, allowing the 
phase of the order parameter to rapidly change by 2ir - 
give rise to residual resistance and can even quench su- 
perconductivity completely. The tiny cross-section of ID 
nanowires reduces the free energy barrier arising from 
a loss of condensation energy for the creation of phase 
slips. Thermal activation of phase slips (TAPS) across 
this barrier is responsible for the residual resistance just 
below T c (LAMH picture) 0. On the other hand, macro- 
scopic quantum tunneling of phase slips (QTPS) as the 
source of residual resistance at low temperatures has re- 
mained controversial despite intense experimental effort 
@; B 0j HI, S 0i S| ■ The observation of macroscopic quan- 
tum tunneling is of significance not only for ID super- 
conductivity, but also for understanding the decoherence 
of quantum systems due to interaction with their envi- 
ronment. Here we clearly establish the presence of QTPS 
in superconducting (SC) Al nanowires, showing the dra- 
matic effect they have on the voltage-current (V-I) char- 
acteristics of the wire. 

Phase slips in ID SC nanowires have traditionally been 
studied by measuring the linear resistance as a function 
of the temperature through the normal-superconducting 
(N-S) transition. Any resistance in excess to that pre- 
dicted by thermal activation is attributed to macroscopic 
tunneling of phase slips @, 0, Q, [H, IE 0]- However, this 
attribution can be flawed because of weak links in the 
wire resulting from inhomogeneities Moreover, the 

fitting of the excess residual resistance to quantum-phasc- 
slip expressions often necessitated an ad hoc 

reduction in the free-energy barrier. Recently Rogachev 



et al. [8|] examined the nonlinear V-I dependence of 
Moo.79Geo.21 SC nanowires: they found that the deduced 
residual resistance, spanning 11 orders in magnitude, fol- 
lowed the prediction of classical TAPS alone, thus con- 
tradicting previous claims of QTPS behavior based on 
residual linear resistance by the same group 0, Q . 

In this Letter, we examine the nonlinear V-I character- 
istics and linear resistance of long aluminum nanowires, 
the narrowest one with dimension 5.2 nm x 6.1 nm x 
100 /im (21 x 24 x 400,000 atoms). Our wires are much 
longer than those of similar cross section reported in the 
literature (~ 0.5^m long) [3, IE HI and the ratio of the 
low temperature SC coherence length to the width (or 
height) is also much larger, thus placing us in a regime 
distinct from previous works. Specifically, we study the 
residual resistance deduced from the nonlinear V-I char- 
acteristics in the superconducting state below the critical 
current. Our main finding is that the V-I dependence 
and residual resistance are inconsistent with the classical 
LAMH behavior, but instead are well described by quan- 
tum expressionSj_either derived from an extension of the 
classical model [2J or in a recently proposed power-law 
form 0, [l(| [H| • The good fits to the different quantum 
expressions, which closely overlap each other, are corrob- 
orated by measured residual linear resistance, achieving 
full consistency within the quantum scenario using a sin- 
gle set of fitting parameters. Our results demonstrate 
the importance of non-classical, quantum phase slips in 
ultranarrow, ID superconducting aluminum wires. 

The nanowires were fabricated by thermally evap- 
orating aluminum onto a narrow, 8 nm-wide, MBE 
(molecular-beam-epitaxy) grown InP ridge, while at once 
linking and partially covering the four-terminal Au/Ti 
measurement pads (Fig. 1(b)). The fabrication is de- 
scribed else where Magnetoresistance at 4.2 K (above 
T c ) was used to characterize the wires (Table 1). The cal- 
culated dirty limit £ [3 (« 94 nm for si, ~ 128 nm for 
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TABLE I: Parameters for sample si and s2, including L 
(length), w (width), t (thickness), Rn (normal state resis- 
tance), p (resistivity), Z e (mean free path), £ (GL coherence 
length), and J c (critical current density), w, t, p and l e were 
deduced based on magnetoresistance measurements [l|. <; 
was calculated using £ = 0.85\/£c>'e where £o = 1600 nm is the 
GL superconducting coherence length in bulk aluminum. 
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FIG. 1: (Color) (a) Linear resistance versus temperature for 
wire s2 in different mag netic fields; the pad series resistance 
has been subtracted |15l|. Inset: The pad series resistance in- 
creases at magnetic fields beyond the critical field of the pad 
regions covered by aluminum (T=0.38, 0.80, 0.99, 1.10 K). 
(b)-(d) Schematic and SEM images of a nominally 8 nm wide 
Al nanowire (similar to samples si and s2): the Al layer does 
not entirely cover the Au/Ti pads which were measured in 
series with the superconducting wire. 



s2) far exceeds the lateral dimensions (~ 5 — 9 nm). 

To investigate the SC behavior, the linear resistance 
through the N-S transition was measured at a current of 
1.6 nA (Fig. 1(a)), while the nonlinear V-I character- 
istics was measured in both the constant-current (Figs. 
2(a)-(b) and 3) and constant-voltage modes below T c jl4j |. 
The linear resistance drops by several orders of magni- 
tude from the normal state resistance (Rn) through T c , 
but remains finite. In Fig. 1(a) we show the linear re- 
sistance versus temperature for wire s2 in different mag- 
netic fields (iJ)|l3- At a given temperature, the resis- 
tance increases with increasing H as superconductivity 
weakens. In the constant-current mode, finite residual 
nonlinear voltage is observed in the V-I curves below the 
critical current jump. Such residual voltage is unobserv- 
able in large wires. In the constant-voltage mode, the V-I 
curves exhibit non-hysteretic voltage steps and S-shaped 
curves down to T/T c ss 0.2, typical of ID superconduc- 
tors d, [IB]. Evidences for wire homogeneity include: 1) 
comparable critical current density in wires si and s2 
2) all V-I traces in the constant-current mode showed a 
single critical current jump. 

To assess the importance of QTPS, we focus on the 
residual resistance and residual voltage below T c for sam- 
ple s2; the wider si will be used for control as discussed 



below. The main datasets consist of the nonlinear V-I 
curves in constant-current mode (Figs. 2(a)-(b) and 3 
(black curves)), and the linear resistance versus T. This 
latter resistance is obtained in two ways: (i) from the 
data in Fig. 1(a) [l5| (Figs. 2(c)-(d) (black curves)), and 
(ii) from the fits to the nonlinear V-I curves according to 
Eqs. (lb) and (3) below (A, Figs. 2(e)-(f)). It will be 
demonstrated through quantitative analysis that this en- 
tire set of data is inconsistent with the TAPS scenario, as 
evidenced by the poor fits to the LAMH expressions (Eqs. 
(la) and (lb)) shown in Figs. 2 and 3 (blue curves). In- 
stead the inclusion of the QTPS contributions (Giordano 
(GIO) model Eqs. (2) and (3)) in the fits is essential in 
providing a fully consistent description based on a single 
set of parameters (see red curves in Figs. 2 and 3). 

Thermal activation of phase slips, important at T < 
T c , is well described by the expressions derived by Langcr 
Ambegaokar, McCumber, Halperin (LAMH) [H. fl7j| : 
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where the quantum resistance R q = irh/2e 2 , attempt fre- 
quency Q = (L/^)(AF/k B T) 1 ^ 2 (h/T GL ), free energy bar- 
rier AF = (8V2/3)( J ff t 2 h /87r)AC, and Ginzburg-Landau 
relaxation time t gl = (tt/8) [h/k B (T-T c )] , I = Aek B T/h. 
L is the wire length, T the temperature, ks the Boltz- 
mann constant, H t h the thermodynamical critical field, 
A the wire cross section, and £ the GL coherence length. 

Quantum tunneling of phase slips, which takes place 
in parallel with thermal activation, is expected to domi- 
nate at low T. GiordanojH proposed a quantum form by 
replacing kT with (1/ a GIO )(h/T GL )^^: 



L AF ( AF 

Rgio = Rq-p \ a-Gio cxp ~a GIO —, 

C V ri/T GL \ n/T GL 



(2) 



where a G iO is a numerical constant of order unit y an d, in 
analogy with the thermal case, we propose that[17l. 1 1 SI] : 



Vqtps=IgioRgio smn (I/Igio) 



(3) 



where I G io = 2e/ {^r GL a GIO )- In this quantum case, a 
resistance similar to the Giordano expression (and nu- 
merically equivalent) has been recently derived, on a mi- 
croscopic basis, by Golubev and Zaikin Q and by Khleb- 
nikov and Pryadko [Io| . In the nonlinear regime, both 
theories predict a crossover from an exponential depen- 
dence to a power- law behavior (V cx I v ) at low T. 

To begin the analysis, we first extract the parameters, 
T c and £ at each magnetic field and a single a G io, by 
fitting to the linear R versus T traces in Figs. 2(c)-(d); 
T c and a G io are input parameters in the analysis of the 
V-I curves. Below T c , a wire is modeled as a normal and 
a superconducting wire in parallel, the latter with a resis- 
tance produced by the relevant phase slips mechanism 
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FIG. 2: (Color) Nonlinear V-I curves and linear resistance 
for sample s2 at different magnetic fields (H): Black curves- 
data; red curves-fits to the GIO (=GIO + LAMH) expres- 
sions for QPS (=TAPS + QTPS), and blue curves-fits to the 
LAMH expressions for TAPS alone, (a)-(b) V-I curves off- 
set for clarity. The fits to QPS are of higher quality com- 
pared to TAPS; each fit includes a series resistance term 
Vs- These V-I curves can be fitted equally well by a power 
law form V = Vs + K ■ (I/Ik)" where 12 > v > 3.2 for 
< H < 1.05 TpjJEl. ( c )-( d ) Linear resistance after back- 
ground subtraction (see Fig. 1(a) and [II]). The LAMH fits 
are poor at low T. (e)-(f) The resistance contribution due 
to phase slips {Rqps) extracted in (a)-(b) from fits to the 
V-I curves using the GIO expressions (discrete points, A). 
Rqps and the linear resistance from (c)-(d) are refitted using 
the GIO expressions (red) with the same aaio(= 1-2) while 
disregarding the irrelevant shoulder feature (dashed line). 



To account for the shoulder feature at a resistance « 1/20 
of i?jv, we assume our wire to be composed of segments of 
two slightly different cross sectional with a fixed 

ratio of 90% (AF oc ^4). The thinner segments corre- 
sponding to the smaller cross sectional area have a larger 
T c ffy- These segments, with critical temperature (resis- 
tance) Tei (i?i), are responsible for > 95% of Rjy. At zero 
field, the fitting parameters are T c i(H = 0), T C 2(H = 0), 
i(H = 0), R u and a GIO , while at H > 0.465 T, only 
T cl (H), T c2 (H), £(H) are varied. Fitting using two £s 
produced unphysical £ dependence on H and was dis- 
carded. The fits to the GIO (red) and LAMH (blue) 
theory are presented in Figs. 2(c)-(d). At H = T both 
fits reproduced the data over several decades in resistance 
and are nearly equivalent above «5 tt. At higher field the 
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FIG. 3: (Color) The nonlinear V-I curves, at T/T c < 0.5 are 
refitted with V = Vqtps + Vs and V = Vtaps + Vs (see 
text).(a)-(c) The data (black), QTPS (red) and TAPS (blue) 
fits and are shown after subtracting the linear background. 
Dotted (dashed) line in 3(a) is the fit to the TAPS expression 
with a 50% reduction (increase) in temperature, respectively. 



GIO theory better models the data, particularly at low 
Tpjj]. The fitting parameters T c \ and £ (o) are presented 
in Fig. 4 (a)-(b). Their H dependence can be fitted to 
simple theoretical expressions as shown. The single value 
of clqio =1.2 we obtain, of order unity as expected, will 
be used throughout all the subsequent analysis. 

Despite the success in fitting to the linear resistance, 
it is not possible to completely rule out weak links as 
the source of residual resistance. To clearly demonstrate 
the importance of quantum phase slips, it is necessary 
to analyze the nonlinear V-I dependence. Generalizing 
the nonlinear analysis of Rogachev et al. Q by including 
QTPS while also accounting for series resistances, the to- 
tal voltage drop across the superconducting nanowire is 

V = Vtaps + Vqtps + Vs, where Vtaps and Vqtps 
are given in Eqs. (lb) and (3), respectively [l7j, and 
Vs = Rsl- Rs is a series resistance which includes the 
contribution of the pads and other ohmic-like contribu- 
tions such as proximity effect of the normal pads on the 
SC wire [2ol . The fitting parameters are Rtaps, Rqtps 
and i?s[17j], while aaio{= 1-2) and T C (H), which enters 
through t gl , arc taken from the previous fits (Fig. 4(a)). 
At all temperatures and fields, the fitting curves repro- 
duce the data extremely well, as shown in red in Figs. 
2(a)-(b) and 3^. At higher T (T/T c > 0.7) the contri- 
bution from QTPS is negligible compared to TAPS. At 
low T {T/T c < 0.5) TAPS is expected to be exponen- 
tially suppressed (Eq. (la)) and QTPS should dominate; 
this expectation is confirmed by the fits which yielded a 
negligible Rtaps compared to Rqtps [H| . 

To further substantiate the importance of QTPS at 
low T while at the same time rule out TAPS as a sig- 
nificant cause of phase slips, we directly compare fits of 
the nonlinear V-I curves using V = Vqtps + Vs and 

V = Vtaps + Vs- The fitting parameters for QTPS are 
Rqtps, Rs, and for TAPS they are Rtaps, Rs- The 
curves are presented in Figs. 3(a)-(c) after subtracting 
the linear background for several H > T. Although the 
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FIG. 4: (Color) GIO(=GIO+LAMH) fitting parameters 
to the linear resistance curves, including those shown in 
Fig. 2(c)-(d), plotted versus H. (a) T cl (o) and T c2 (not 
shown) fit well to the pair breaking perturbation theory [TJ 
(green), (b) £ (o) is fitted to both£ = £ c /,yi - (H/H c ) 2 (ma- 
genta) (Ace ^/l - {H/H c ) 2 and £ocl/\/A[I3), and to an ad 
hoc expression £ = £ c /^/l - (H/H c ) 2 (green); £ c = £(// = 0) 
and He is the critical field. £ is better described by the ad hoc 
expression. £ (A) derived from the fits to Rqps, including 
those in Fig. 2. (e)-(f), is fitted to theory as before (dashed 
lines) yielding £ c ~ 135 nm, within 5% of the value in Table 1. 



two expressions have the same number of parameters, the 
QTPS fits (red) are of good quality while the TAPS fits 
(blue) are evidently much poorer. Varying T by ±50% in 
the TAPS fits still failed to reproduce the data, as shown 
in Fig. 3(a) for T = 0.59 K and H = 0.465 T. Comparing 
the x-squares of the fits that include QTPS versus TAPS 
alone using the statistical F-test yielded strong support 
for the importance of quantum phase slips with a confi- 
dence level >99.99%(23. 

The Rqps = Rqtps^Rtaps extracted from the non- 
linear V-I fits enables to exclude contributions to residual 
resistance that arc irrelevant to phase slips, providing a 
means to check full consistency. In Figs. 2(e)-(f) we plot 
Rqps (A symbols) for T <T C and replot the linear resis- 
tance through the N-S transition for comparison; we then 
refit the linear resistance expressions (GIO+LAMH, Eqs. 
(la) + (2)) to these data, keeping aaio = 1-2 unaltered 
while varying T c and £ . The discrete Rqps data points 
versus T no longer exhibit a shoulder feature, and only 
one T c was needed. Satisfactory agreement is achieved as 
depicted in Figs. 2(e)-(f) (red). Similar agreement com- 
parable to the H — 0.465 T trace is obtained at higher 
H = 0.7 T, 0.93 T, etc. T c was essentially unchanged 
compared to T& in Fig. 4(a) but £ (Fig. 4(b), A) was 
reduced at higher H, compared to the data obtained pre- 
viously (Fig. 4(b), o). Fitting to the new values yielded 
£ c = £(H = 0) w 135 nm, within 5% of the calculated 
value in Table 1. As a final consistency check, analysis of 
the wider superconducting sample sl[12| did not produce 
any sign of QTPS for H J$ 0.9 T. This is sensible because 
of the larger free-energy barrier. 

The consistency achieved in our analysis of the nonlin- 
ear V-I curves and residual resistance, supported by the 
quantitative statistical F-test, demonstrates the impor- 
tance of quantum phase slips. This helps rule out other 



scenarios, e.g. TAPS alone or weak links. Furthermore, 
our results establish that the transport properties of ID 
superconducting nanowires at temperatures much below 
T c are determined primarily by the macroscopic quantum 
tunneling of phase slips. Our findings pave the way to 
the study of newly predicted quantum phase transitions 
in metallic nanowires 24 1 . 

F.A. thanks M.E. Rizza for her support. Work sup- 
ported by NSF DMR-0135931 and DMR-0401648. 
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Erratum: Evidence for Macroscopic Quantum Tunneling of Phase Slips in Long 
One-Dimensional Superconducting Al Wires 
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1 Physics Department, Duke University, Durham, NC 2770^\ 
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(Dated: February 2, 2008) 

PACS numbers: 74.78.Na, 74.25.Fy, 74.25. Ha, 74.40.+k 

The definition of t gl should read t gl = (ir/8)[h/k B (T c - T)]. 

An inadvertent error in fitting to the magneto-resistance of the Al wires led to an underestimate of the wire widths 
and cross sections Q. The revised width/cross section for the Al wires si (10 fj,m long) and s2 (100 /im long) are 
respectively: 11.4 nm /130 nm 2 and 7.5 nm/58 nm 2 roughly an increase of 40 % in width and factor of 2 in cross 
section. The wire widths have an error of ± 10 %, and the cross section error is ± 20 %. The revised sample 
parameters are listed in the revised Table I. Note that the wire length, L. and normal state resistance, Rn, remain 
unchanged. 

The dependence of the free-energy barrier AF on the cross sectional area A is contained in the quantity i/i?Ar0], 
removing any explicit dependence on A. Thus, the revised quantities, including A, are not used in the main analysis, 
and do not affect the main conclusion of the paper. The only notable effect is that the agreement between the fitted 
coherence length £ c and the calculated £ now occurs at the boundary of the error bars. Note that the error on £ c is 
15%. The fitted values for AF{H = 0, T = 0) and £ C (H = 0,T = 0) (135 nm), combined with the revised A, yield a 
thermodynamic critical field H th (H = 0,T = 0) w 0.0096 T (96 Oersted), close to the known value of 0.0105 T (105 
Oersted) for bulk aluminum. 
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TABLE I: Parameters for sample si and s2. The parameters L (length), w (width), t (thickness), Rn (normal state resistance), 
p (resistivity), l e (mean free path), £ (GL coherence length), and J c (critical current density), for the two samples si and s2. 
w, t, p and l e have been deduced following the analysis in Altomare et al.[l| based on magnetoresistance measurements. £ was 
calculated using £ = 0.85\/£o'e where £o = 1600 nm is the GL superconducting coherence length. 



Evidence for Macroscopic Quantum Tunneling of Phase Slips in Long 
One-Dimensional Superconducting Al Wires 

Fabio Altomare, Albert M. Chang, Michael R. Melloch, Yuguang Hong, and Charles W. Tu 



1. POWER LAW 



An alternative fit to the V-I curves, distinct from the GIO(= GIO+LAMH) fit (see text and Figs. 2(a)-(b) of main 
manuscript), is performed using the power law dependence proposed for quantum phase slips [1]: 



where Vs — Rsl is due to the series resistance Rs, and K, Ik are appropriate coefficients. These fits are indistin- 
guishable from the best fits with the GIO model and will not be reproduced here. The exponent v is plotted as a 
function of the magnetic field H in Fig. El, and is obtained from the V-I curves at different temperatures in Figs. 
2(a)-(b) and 3 (main manuscript) and from additional H that are not shown. The large exponent value (> 3.2) and 
the rapid decrease of v with increasing H is consistent with naive estimates based on the theoretical model. 

Furthermore, it is also possible to account for the observed linear resistance versus T behavior with a power law 
dependence on the temperature: Rkl = K\T V . Theory [1, 2] predicts that the power-law exponent v in the R 
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FIG. El: The exponent v, extracted from the fit of the V-I curves to the power law expression, has been plotted as function 
of the magnetic field. Different symbols correspond to V-I curves at different temperatures. 
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versus T dependence and in the V-I curves should be identical at each magnetic field. However, K\ is predicted only 
under special limiting conditions. The fitting therefore necessitates the addition of K\ as a fitting parameter. In 
doing so, the fitting becomes less restrictive. 

2. DETAILS ON THE SUBTRACTION OF THE SERIES RESISTANCE 

In the raw data for the 4-terminal, linear resistance versus temperature trace for sample s2 measured at H = T, 
the resistance decreases precipitously from its normal state value, Rn = 86fcSl (3> R q = hj '4e2 sa 6.45fcSl), through the 
superconducting transition and saturates below 0.75 K at 47± 4 SI. This saturated value of the resistance is due to 
the pads and is measured in series with the nanowire. There are two contributions to the series resistance from each 
pad: one from an Au/Ti region without aluminum (Al) coverage which remains after the final aluminum evaporation 
(47 SI), and the other from the covered Al/Au/Ti region (0 SI when superconducting, and « 423 SI when normal). 
The superconducting-to-normal (S-N) transition of the Al/Au/Ti pad regions as a function of the magnetic field H is 
shown in the inset of Fig. 1(a) in the main manuscript for several temperatures. It is indicated by the resistive jump at 
H ~ ±0.1 T. To confirm that the pad is responsible for this jump, separate two terminal measurement was performed 
for each pad, and in addition a different sample with the same Al/Au/Ti trilayer pad geometry, but connected to a 
wide Al wire, was also measured. In both cases a similar jump was observed in the resistance near H ~ ±0.1 T. At 
0.34 K, the total pad resistance becomes independent of H above H m 0.35 T, having attained the full normal state 
value of w 470 SI. Above this field, this normal state value is also temperature independent . To extract the resistance 
of the nanowire, we subtract this resistance of 470 SI from the R vs. T raw data traces when H > 0.465 T, while for 
H < 0.03 T we subtract 47 SI instead, as the Al/Au/Ti regions become superconducting below «1 K and no longer 
contribute to the series resistance. The resultant traces are plotted in Fig. 1(a) of the main manuscript. 

3. F-TEST 

In the nonlinear voltage due to the thermal activation of phase slips (TAPS), Vtaps = IoRLAMHsinh{I / Iq) 
(Eq. (lb) main manuscript), the current scale Iq = AeksT/h does not depend on the magnetic field H. However in 
the nonlinear voltage due to quantum tunneling of phase slips (QTPS), Vqtps = lGioRGio s i n h{I / Igio) (Eq. (3) 
main manuscript), the current scale Igio = ^/(^tglo-gio) does depend on H through t gl = (tt /8)[H/ ks{T — T C (H))]. 
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Both T C (H) and acio (in Igio) are deduced from fits to the linear R versus T traces and their values are presented 
in Fig. 4(a) and the text of the main manuscript, respectively. 

In the measured voltage-current (V-I) curves the importance of TAPS was tested by fitting the experimental curves 
to V = Vtaps + Vs- The series voltage Vs = Rsl represents the contribution from ohmic-like resistances in scries 
with the superconducting wire. To account for the contribution from the quantum tunneling of phase slips (QTPS), 
a Vqtps term is added. The total QPS (= QTPS + TAPS) contribution is given by V = Vqtps + Vtaps + V s . 
The fitting parameters for the TAPS contribution without QTPS are Rtaps and R$ [3]. For the QPS fit including 
QTPS, a third parameter Rqtps is added. Note that to distinguish the fitting parameters from the theoretical ones, 
we refer to Rtaps (Rqtps) as the TAPS (QTPS) resistance extracted from the fits to the V-I curves, and to Rlamh 
(Rgio) as the theoretically calculated values (see Eqs. 1(a), 1(b), 2, and 3 in the main manuscript). As can be seen in 
Fig. 2(b) (main manuscript), the QPS fits are of good quality while the TAPS fits are much poorer. To quantitatively 
compare the two different fitting models, the statistical F-test[4] on chi-squares was applied to the random variable 
F(n, m) = [(xqps ~ Xtaps) I X Q ps\l[(dfQPS ~ df T APs)/df Q Ps\. Here n = df QPS - df T APS, m = df QPS with df the 
number of degree of freedom, in the model under exam, defined by df = N — p, where N is the number of points 
and p the number of fitting parameter in the model, and Xqps(Xtaps) represents the chi-square obtained by fitting 
the experimental data to the QPS (TAPS) model: x 2 (ai---a n ) = -^^[yi — f(xi, d\ ■ ■ ■ a n )], with j/j and X{ the 
experimental data and a\ - ■ ■ a n fitting parameters. The F-test is useful for comparing two fitting expressions which 
belong to the same family of related expressions but one with more fitting parametcr(s) than the other. 

To accept the QPS and refuse the TAPS model with a confidence level exceeding 99%, F(n,m) has to exceed the 
criterion for an integrated probability of 99% in the F(n, m) distribution with n and m degrees of freedom. The F-test 
indicates that for the low temperature traces (T/T c < 0.5) displayed in Fig. 2(b) (main manuscript) as well as for all 
experimental V-I curves in Figs. 3 (main manuscript), the QPS model is prcferred[5] with a confidence level greater 
than 99.99%. For example, for the traces in Fig. E2 (from Fig. 3(a) of the main manuscript), with: n=l, m=89, 
Xtaps = 1-13-E 1 — Hi Xqps = 1-82-E — 12, we obtain that the likelihood in favor of the QPS model is basically 100% 
(actual numerical value 1 — 2.8 x 10~ 13 ). Moreover, the excellent QTPS fits (no TAPS term) presented in Fig.3(a)-(c) 
(main manuscript) indicate that TAPS does not contribute in the low temperature limit (T/T c < 0.5). 
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FIG. E2: (Color) Even by increasing the temperature by as much as 50% the TAPS model does not properly reproduce the 
experimental data as well as the QPS model (shown here for T=0.59 K and i/=0.465 T). Reducing the temperature makes the 
fit even worse. The traces are displayed after subtracting the linear background contribution due to the series resistance. 



In this section of the EPAPS we will show that any heating effect that may arise from the non-negligible power level 
dissipated in the superconducting wire, has minimal effects in the analysis of the nonlinear V-I curves of Fig. 2(a)-(b) 
and 3 (main manuscript) below the critical current Ic (I < Ic) an d can be neglected, without loss of generality, 
in our fitting procedure. We will also show that it is reasonable to attribute the dominant thermal conduction to 
the lattice conduction of the InP ridge. After providing some preliminary details, we will show that, at a power 
of 15 pW dissipated through the nanowire and at a bath temperature Tb at ft=0.4 K, this self-heating would increase 
the wire temperature (T w ) to no more than 0.48 K (a relative change in temperature AT/T < 20%). Additionally, 
this temperature rise decreases with increasing bath temperature, due to increased thermal conduction of the InP 
ridge on which the nanowire sits, and, under the same assumptions, at T=0.6 K the relative change in temperature 
would be AT/T < 4% and at T= 0.7 K AT/T < 2.3%: this change in temperature alone, which does not depend on 
the magnetic field, cannot account for the observed non-linearity in the V-I. Our situation is in direct contrast with 
the results presented in other works where the superconducting wire was fabricated on suspended nanotubes[6, 7] 
and therefore lacked a direct thermal link to the substrate aside from the wire itself, for which, the electron thermal 
conductivity deduced from the electrical conductivity via the Wicdemann-Franz relation, was the dominant (but 
inadequate) means for thermal conduction. In particular, in the work of Johannson et al., heating seriously limits 
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the range of viable currents at their measuring temperature (~ 10 mK, much lower than in the present work). The 
two systems, however, cannot be directly compared because in the suspended nanotube geometry, the only additional 
cooling effect, beyond the electric thermal conduction, is due to the 3 He bath which is negligible because of the Kapitza 
boundary resistance. In our case the InP ridge, on which the wire is fabricated, readily thermalizes with the wafer; this 
in turn, because of its large (compared to the nanowire and the ridge) physical size (0.5 mm x 1.5 mm x 2 mm) is not 
affected by the Kapitza boundary resistance and can be efficiently cooled by the bath. We will demonstrate that, a) 
assuming a constant effective wire temperature (T w ), but higher than the bath temperature (Ti, a th), the TAPS theory 
still cannot properly reproduce the experimental data; b) the small increase in temperature we estimate, and the 
associated increase in resistance, is not nearly sufficient to explain the observed non linearity and, even by modifying 
our fitting procedure so as to include this temperature change into a temperature dependent coefficient, the TAPS 
model does not adequately describe the experimental data. On the other hand, the QPS (= QTPS + TAPS) model 
fits the data well at all temperatures and fields examined. We would like to point out that at T/T c < 0.5, which will 
be the main focus of this section, in the QPS model Vtaps is negligible with respect to Vqtps- Therefore, in the 
following we will refer equivalcntly to "QPS model" or to "QTPS model" . 

The first point (a) is readily proven in Fig. E2 where the experimental data (for a representative trace displayed in 
Fig. 3(a) of the main paper) are fitted with the TAPS theory by changing the temperature by as much as 50%: the 
agreement between the TAPS model and the experimental data is still poor when compared to the QPS model. 
Regarding the second point (b), we make the following observations to render the discussion more quantitative: 
1) The power dissipated in the wire can be deduced from Fig. 2(a)-(b) using V ■ I = R ■ I 2 after subtracting the 
contribution due to the pads — 50 fl at H=0 T and 470 f2 at H > 0.465 T. Immediately before switching from the 
superconducting state to the normal state, the power dissipated in the wire (at all fields considered) is in the range 
P f» 5-15 pW. To see that the electronic thermal contribution alone is not sufficient to cool the wire, as pointed out by 
referee A, we use a simple estimate based on the Wiedemann- Franz law and the heat diffusion equation (Q = An\7T) , 
with A the cross sectional area of the wire, k the thermal conductivity and VT the temperature gradient, to find 
that, at T=0.5 K, it would produce an increase in temperature of several kelvin: an entirely inaccurate scenario, since 
the nanowire still exhibits a non-zero critical current, and therefore the effective wire temperature is smaller than 
the critical temperature (T c ~ 1.7 K). This can be seen most dramatically in the V-I curve in the down sweep of 
the current (Fig. E3, red curve), where the power in the wire, immediately before switching from the normal to the 
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FIG. E3: (Color) V-I characteristic for the up (black) and down (red) sweep at T=0.6 K and _H=0.465 T. The dashed line 
intersects the current axis at a finite current indicating the presence of a nonzero average supercurrent. 

superconducting state, is greater than 200 pW. If there were no additional channels, i.e. the InP ridge and substrate, 
to conduct away the heat generated by this large power, the temperature rise would be so huge (> 10 K) that no 
transition from the normal to the superconducting state would be observed. In fact, based on the experimental data, 
even at these power levels, the wire temperature cannot be above the critical temperature because a nonzero average 
supercurrent still flows in the nanowire. This can be deduced by noticing that the extrapolation of the down sweep 
V-I intercepts the current axis at a finite value (a typical curve is shown in Fig. E3 for the V-I trace taken at T=0.6 K 
and ff=0.465 T). While this does not exclude the possibility of a slight temperature rise, it clearly excludes a naive 
model of thermal conduction with dominant electron component, based on the Wiedemann- Franz law, which would 
yield a large change in temperature. Moreover, this failure implies that the above description does not adequately 
model the experimental situation. 

2) Thermalization of the wire in the vertical direction (~ 6.1 nm), as opposed to along the wire through the pads 
can readily be achieved via electronic conduction. This means that heat can easily be transferred from the top to 
the bottom of the nanowire across the 6.1 nm distance (the farthest and closest to the ridge, respectively), and any 
difference in temperature quickly equilibrated. In fact, because of the large cross sectional — 100 nm x 5.25 nm — 
surface area (versus 5.25 nm x 6.1 nm cross section for thermal conduction along the nanowire), even for a very 
high power of P~ 300 pW (which would render the wire normal) and T=0.5 K, by using the heat diffusion equation 
combined with the Wiedemann-Franz law, we obtain AT < 1 • 10 -5 K. Any heat generated in the wire, therefore, will 
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be quickly transferred to the InP ridge which (as it will be discussed later in detail) is the dominant mechanism for 
thermally anchoring the wire to the bath. 

3) The crystal thermal conductivity (k) of the InP ridge and the InP/GaAs substrate is independent of the magnetic 
field and k oc C v X p h where X p h is the phonon mean free path and C v oc T 3 is the specific heat. For the bulk InP/GaAs 
substrate the phonon mean free path is limited by boundary scattering while for the InP ridge is most likely cut off by 
the height of the ridge (~ 25 nm). Since the ridge width is only ~ 8 nm and is almost atomically flat, the reflection off 
the ridge walls is likely specular (note that the thermal wavelength at 1 K is sa 50 nm). Only bulk phonons, however, 
will be effective in transferring the heat because surface defects and adsorbates will limit the heat transferred by 
surface phonons. 

4) One of the biggest problems in thermally anchoring a sample at the bath temperature, in low temperature 
measurements, is the Kapitza boundary resistance, defined as 

R k = AAT/Q (2) 

where A is the contact area, AT the difference in temperature between the bath and the sample (nanowire in our 
case) and Q the heat flow. Based on values for other typical metals [8], because of the small contact area between 
3 Hc and the aluminum nanowire, the bath cannot properly cool the nanowire. On the other hand, the crystal size 
(.5 mm x 1.5 mm x 2 mm) is such that any thermal boundary resistance can be ignored. Additionally, given the 
typical lattice conductivity of good crystals (at T=l K, k r* 0.01 — 8 W/K m)[9, 10] from the heat diffusion equation, 
at T=10 pW, the increase in temperature for the crystal would be AT < 1 ■ 10~ 4 K. Based on the above discussion, 
since the aluminum nanowire exhibits superconducting behavior, even at power of the order of 200 pW or more (finite 
non-zero switching current), the temperature must still be below T c ~ 1.7 K. Therefore there must be some channel 
that efficiently anchors the nanowire to the bath temperature: this channel is the InP ridge, directly connected to the 
crystal substrate. 

Below we present a quantitative analysis of the temperature rise, at the relevant power levels of P < 15 pW, for the 
non-linear V-I curves for Figs. 2(b) and 3 of the main manuscript. We begin by noting that the change in resistance 
due to the heat is not sufficient in explaining the observed non linearity. 

In fact, let us assume (this will be demonstrated later) that a power of 15 pW produces a relative increase in 
temperature of 20% at T=0.4 K, 4% at T=0.6 K and 2.3% at T=0.7 K. Because of the temperature dependence of 
the resistance-i. e. it increases with increasing temperature as shown in Fig. 1(a) of the main paper, there will be 
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FIG. E4: (Color) Experimental V-I curves at _ff=1.16 T (black) and voltage drop across the superconducting nanowire calculated 
(red) up to the switching current by attributing the resistance change to the change in temperature (from bottom to top: 
AT/T = 20% at T=0.4 K, AT/T = 4% at T=0.6 K, AT /T = 2.3% at T=0.7 K). The contribution due to the series resistance 
linear in the current has been subtracted and o indicates curves multiplied by 4. Even after increasing our estimate by an 
additional 10% (green) the agreement is still poor. Curves are offset for clarity. 

nonlincarity in the V-I characteristics. Any heating effect, however, will be independent of the magnetic field, so we 
can use the R vs. T (in the linear regime, Fig. 1(a) main manuscript) curve at £f=1.16 T as a thermometer. This has 
the advantage of providing, in the temperature range of interest, the biggest change in resistance as function of the 
change in temperature. Since at this magnetic field the R vs. T data stops at T ~ 0.48 K, we linearly interpolated 
to obtain the data at the lowest temperature (T=0.42 K) by using the resistance in the corresponding V-I curve (at 
H=1.16 T) deduced in the limit of low current. Our results are shown in Fig. E4, where both experimental data 
and the voltage drop across the superconducting wire, calculated from the measured resistance with the change in 
temperature taken into account, are displayed after subtracting the contribution due to the series resistance which 
is linear in the current. Increasing the temperature by an additional 10% still does not reproduce the experimental 
data. 

Therefore we can conclude that the major portion of the non-linearity cannot be attributed to self-heating. 

From our discussion above, it is clear that assuming the wire temperature either does not change or changes within 
our estimate, the QPS model better describes the experimental data (see also Fig. E9). The only remaining task 
is to describe the method for obtaining the previously stated upper-bound estimates (AT/T = 20% at T=0.4 K, 
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FIG. E5: (Color) //=0.465 T: (a) the down sweep at T=0.7 K has the same switching current as the up sweep at T=1.04 K: we 
conclude that at Tb at h~0.7 K, a power level of P ~ 600 pW raises the wire temperature up to 1^ = 1.04 K. (b) At T=0.39 K, 
the down switching current is in between the up switching current at T=1.04 K and T=0.89 K. Linearly interpolating between 
these two values (see also Fig. E6), we obtain that at Tbath ~ 0.39 K, a power level P ~ 644 pW increases the wire temperature 
to T w = 1.015 K. The dotted curves refer to up and down sweep not relevant to this discussion. Inset: the full range of the V-I 
curves for the three different temperatures is shown. 

AT/T = 4% at T=0.6 K, AT/T = 2.3% at T=0.7 K), which will involve using the position of the critical current 
jumps in the V-I curves themselves. Since self-heating will be particularly important at lower temperature, we will 
concentrate on the curves presented in Fig. 3 of the main manuscript for which the difference between the QPS (or 
the numerically equivalent QTPS) model and the TAPS model is the most apparent. Let us consider the situation 
depicted in Fig. E5(a). Here we show for magnetic field H = 0.465 T up (black) and down (red) sweeps at a bath 
temperature of T=0.7 K together with up (green) and down (blue) sweeps at T=1.04 K. It is evident that the current 
at which the nanowire switches from the normal to the superconducting state in a down-sweep at the bath temperature 
Tbath=0-7 K (red), with a corresponding dissipated power of P ~ 600 pW, coincides with the current at which the 
nanowire switches from the superconducting to the normal state at a bath temperature Tb at ^=1.04 K in an up-sweep 
(green). Since the switching positions overlap, we infer that at a bath temperature Tb a th~0-7 K the wire temperature 
is not more than 7^= 1.04 K, despite a power level P ~ 600 pW. This is equivalent to saying that, with a power 
P ~ 600 pW dissipated in the nanowire and at a bath temperature 21,^=0.7 K, the wire temperature rises to no 
more than T W =1.QA K. This identification can be made because in the up-sweep, immediately before the transition 
from superconducting to normal state, the power level is sufficiently low that the wire temperature is very close to 
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FIG. E6: At a field H=0A65 T, the bath temperature has been recorded as function of the up (black) and down (red) switching 
current. If we were to perform a measurement at T=1.015 K, the up switching current would be the same as the down switching 
current at T=0.39 K at a power given by the product V ■ I immediately before the transition from normal to superconducting 
state. By repeating the same procedure for the other fields of interest (H=0, 0.7,0.93 T), we obtain the datapoints in Fig. E8 
after allowing for some uncertainty in the value of the down switching current. 

the bath temperature. The 600 pW dissipated power is estimated by multiplying the current and the voltage at 
the switching current in the down sweep (at Tf, at ^=0.7 K) after subtracting the contribution due to the pads. For 
situations where a good match between different bath temperatures is not available, we linearly interpolate between 
two nearby temperatures to extract the value of interest (Fig. E6). An example is provided in Figs. E5(b) and E6. 
At a bath temperature Tb a th—0.39 K, the wire temperature is T„,=1.015 K at a power level P ~ 644 pW (as before 
the power level is calculated from Fig.E5(b) after subtracting the pad contribution). 

By following this procedure for all our datasets (representative VT are shown in Fig. E7), we obtain, for every T^ath 
of interest (T=0.4 K, 0.6 K and 0.7 K), the corresponding T w as a function of the power dissipated in the nanowirc. 
This allows us to estimate, for the lattice temperatures mentioned above, the temperature rise due to the dissipated 
power: in particular we will be interested in the change in temperature at the power levels at which the wire switches 
from the superconducting to the normal state (< 15 pW). 

Our results are displayed in Fig, E8 (a)-(c) where the scatter in the datapoints at a given magnetic field reflects 
the uncertainty in the exact determination of the switching current. In order to analyze these data, we model the 
InP ridge thermal conductivity as a 1-dimensional problem with P = Q = kAVT. By integrating the heat diffusion 
equation along the ridge height direction and taking into account the the temperature dependence (T 3 ) of the crystal 
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FIG. E7: (Color) (a)-(d): Typical V-I characteristics up (black) and down (red) sweep at T « 0.6 K at different magnetic field. 



thermal conductivity, in a scenario where \ p h ~ const, the expected behavior of T w is T w = P/C + T? ' . , with C 
an appropriate constant that depends on the crystal and on the ridge geometry. Note that at zero power dissipated 
through the nanowire (and therefore zero current flowing), the aluminum nanowire must be in thermal equilibrium 
with the bath. 

From the fit to this functional form, displayed in Fig. E8, we find that, for P=15 pW, AT/T = 20% at T=0.4 K, 
AT/T = 4% at T=0.6 K and AT/T = 2.3% at T=0.7 K. 

It must be noted that these values are an upper bound because we have neglected the temperature dependence 
of the phonon mean free path (which may behave like X p h oc 1/T[9]): the thermal conductivity, therefore, might 



have a weaker dependence on the temperature. By assuming T w = y P/C + T^ ath , with a a fitting parameter, we 
obtain a = 2.8 in good agreement with our expectations (a ~ 3), which would further reduce our upper bound 
estimate for AT/T. By incorporating the theoretical temperature dependence into Rtaps and Rqtps based on 
Eqs. (la) and (2) of the main manuscript, and repeating the same analysis that has led to Fig. 3 (main manuscript), 
we obtained the fits shown in Fig. E9. It is apparent that the experimental curves can be adequately fitted only by 
the QPS (or numerically equivalent QTPS) expression. The quality of these fits is comparable to quality of those 
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FIG. E8: (Color) Temperature reached by the superconducting nanowire at different power level and fit with the functional 
form discussed in the text for (a) Tb a th~0A K, (b) Tb a th~0.6 K, (c) Tb a th=0-7 K: the scatter in the datapoints is due to the 
uncertainty in estimating the value of the current at which the nanowire switches from the normal to the superconducting 
state. Different symbol shapes indicate data points obtained at different fields. 

shown in Fig. 3 of the main paper at T. w = Tj, t/,. In Fig. E9 when fitting using the full expression for the QPS model 
(V = Vqtps + Vtaps + Vs)i the contribution from Vtaps ls found to be negligible compared to Vqtps, and F-test 
on the chi-squares for the datasets in Fig. E9 supports the QPS model with a confidence level exceedingly 99.99%. 

As an additional independent check on these temperature estimates, we point out that since the temperature 
increase is due solely to the power dissipated in the nanowire and should not depend on the applied magnetic field, 
we can use the V-I data at 7J=1.16 T to calculate a resistance (V/J) and a corresponding power level at each /, and 
then use the R vs. T curve to convert it to an effective, rough upper-bound temperature assuming that any increase 
in resistance comes entirely from the temperature rise alone and not from other sources of nonlincarity. By doing so 
we would obtain, for a power P=15 pW, AT/T = 31%, 17%, 17% at T=0A, 0.6, 0.7 K, respectively. By taking into 
account this estimated rise in temperature, larger than the more stringent estimate provided above, the TAPS theory 
still would not reproduce the data in Figs. 2(a)-(b) and 3 (main manuscript). 
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FIG. E9: (Color) The curves displayed in Fig. 3 (main paper) have been fit after including the theoretical temperature 
dependence in the fitting parameters Rtaps and Rqtps- As before, only the QPS (or numerically equivalent QTPS) model 
adequately fits the experimental data at all fields and temperatures. (The linear contribution due to the series resistance has 
been subtracted and the curves have been offset for clarity). 

The last issue requiring examination is whether there is a possibility that the electron and phonons can become 
decoupled in the nanowire, which would prevent the transfer of heat from the electrons to the lattice, resulting in 
a situation where the electrons and lattice are no longer in equilibrium. The clcctron-phonon coupling follows a T 5 
dependence in temperature, and is generally believed to be sufficiently strong at T > IK. On the other hand, because 
of the T 5 dependence, a decoupling may arise below T ~ 1 K. We perform the following estimate to demonstrate 
that the reduced coupling at lower T does not cause a problem for heat transfer at power levels < 15 pW. To set 
the scale of the coupling we use the behavior at an electron temperature of 1 K for reference, where the electrons 
and phonons are expected to be in thermal equilibrium at high dissipated power. For a bath temperature of ~ 0.4 
K and electron temperature of 1 K, a dissipated power of 633 pW in the nanowire can be transferred to the lattice 
without difficulty (Fig. E8(a)). At the lower electron temperature of 0.48 K corresponding to our estimated upper 
bound of a 20% temperature rise at 0.4 K bath temperature, the electron-phonon coupling is reduced by a factor 
(T/IK) 5 = 0.0255. The power level which would still allow the electrons and lattice to be in thermal equilibrium 
would thus be 633 pW*0. 0255=16. 14 pW. The power relevant to the nonlinear V-I curves in our analysis (< 15 pW) 
is below this value. Therefore the reduction in coupling should not be problematic. This issue is even less relevant for 
the other two temperatures of interest presented in Figs. E8(b) and E8(c) (bath temperatures > 0.4 K and electron 
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temperatures > 0.48 K). 

In view of the evidence (both theoretical and experimental) , we conclude that while slight increases in temperature 
cannot be ruled out for the power levels under consideration, the increases are inconsequential for the main results of 
the Letter and can be neglected. Thus, in fact we are observing a nonlinearity in the VT characteristics which cannot 
be explained by the TAPS model, and which is instead compatible with the phenomenological model of QPS. 
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